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Abstract 

In this paper we demonstrate that the inevitable action of the environment can be substan- 
tially weakened when considering appropriate nonstationary quantum systems. Beyond protect- 
ing quantum states against decoherence, an oscillating frequency can be engineered to make the 
system-reservoir coupling almost negligible. Differently from the program for engineering reservoir 
and similarly to the schemes for dynamical decoupling of open quantum systems, our technique 
does not require a previous knowledge of the state to be protected. However, differently from the 
previously-reported schemes for dynamical decoupling, our technique does not rely on the availabil- 
ity of tailored external pulses acting faster than the shortest time scale accessible to the reservoir 
degree of freedom. 



1 



A great deal of attention has recently been devoted to quantum information theory owing 
to its strategic position, joining up several areas of theoretical and experimental physics. 
As eventually all domains of low energy physics may provide potential platforms for the 
implementation of quantum logic operations, efforts have been concentrated on overcoming 
some sensitive problems that constitute a spectacular barrier against their realization. These 
problem areas touch on both fundamental physics phenomena — such as decoherence and 
nonlocality — and outstanding technological issues such as individual addressing of quantum 
systems, separated by only a few /im, with insignificant error [l| . 

As the debate around nonlocality seems to be subsiding through a set of experimental 
results — such as i) technological evidence against the so-called loopholes ^, ii) the demon- 
strated violation of Bell's inequality with two-photon fringe visibilities in excess of 97% j^, 
and iii) highly successful experimental quantum teleportation [4| — the program for quan- 
tum state protection is still at an early stage, despite all the achievements. A promising 
suggestion on this subject refers to the possibility of manipulating the system-reservoir cou- 
pling through an additional interaction between the system and a classical ancilla. This 
control of decoherence through engineered reservoirs has beeri^ theoretically implemented 
for atomic two- level systems, exploiting a structured reservoir [5] or mimicking a squeezec- 
bath interaction |6j]. In the domain of trapped ions, beyond a theoretical proposition j?], 
engineered reservoirs have also been experimentally implemented for superposed motional 
states of a single trapped atom j^. Another strategy, also experimentally investigated 
91], involves collective decoherence, where a composite system interacting with a common 
reservoir [l^ exhibits a decoherence- free subspace (DFS). Whereas a common reservoir is 
crucial for shielding quantum coherence in a DFS, the quantum-error correction codes QECC 



ll| work, instead, on the assumption that the decoherence process acts independently on 



each of the quantum systems encoding a qubit. The issue of the physical grounds for the 
assumptions behind a common or distinct reservoirs is in detail in Ref. 121]. 

We also mention a recent proposal for the control of coherence of a two-level quantum 

n n 

system p^, based on random dynamical decoupling methods p^. These methods resemble a 
previous technique to suppress decoherence that used a tailored external driving force acting 
as pulses [isl which, as in the present paper, was applied to a cavity-mode superposition 



state. In Ref. 



16|, in a more general scope, the authors formulated a model for decoupling a 



generic open quantum system from the environmental influence also bailing out on tailored 
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external pulses to induce motions into the system which are faster than the shortest time 
scale accessible to the reservoir degree of freedom. 
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In the present work we achieve the goal of Refs. which goes beyond the 

quest for quantum state protection through engineered reservoir, from a different approach: 
We demonstrate, arguing from quite general and current assumptions, that a nonstation- 
ary resonator could be almost completely decoupled from the environment, rendering the 
damping factor that characterizes the environment negligible. Note that, differently from 
our proposal as well as those in Refs. Q, Q, H, Il6l ]. the schemes of engineered reservoirs 



require a previous knowledge of the state to be protected. Evidently, this requirement for- 
bids the use of engineered reservoirs for the implementation of logic operations, making the 
schemes of switching off the system-reservoir interaction more attractive. Finally, we observe 
that the control of decoherence through the frequency modulation of the system-heat-bath 

n nn 

coupling has been proposed earlier [17|, but as in Refs. [13|, [IJ], such control is achieved for 



a two-level system instead of a cavity mode. 

Assuming a nonstationary mode coupled to the environment, we get the Hamiltonian 

H{t) = u{t)a^a + J2^kblbk + J^^^it) (abl + a^bk^ , (1) 

k k 

with (a) and (6) standing for the creation (annihilation) operators of the nonstationary 
field uj{t) and the kth bath mode Uk, respectively. Assuming the time-dependent (TD) 
relation uj{t) = ujq—x sin(Ct), the system-reservoir couplings also turn out to be TD functions 
\k(t). The simple TD form of the free Hamiltonian Hq = uj(t)a'^a + J^k^^blbk enables us to 
describe, through the transformation Uit) = exp (^—i j^HQ{T)dT^, the Hamiltonian in the 
interaction picture 

V{t) = aK\t)+a^K{t), (2) 

where we have defined the TD operator A(t) = Afc(i(:)6fc exp [iAjt(t)] and parameter 
Ak{t) = ^l(t) — tUkt, with Q{t) = Jl^to{T)dT. For the case of weak system-reservoir coupling 
the evolution of the density matrix of the nonstationary field, in the interaction picture and 
to the second order of perturbation, is given by 

^ = - /^^'Tr/?. [Vit), [Vit'),pniO) ® pit)]] , (3) 

where we have employed the usual approximation Pr{0) ® p{t). Assuming that the reservoir 
frequencies are very closely spaced, with spectral density ^{p), to allow the continuum 



summation of the coupling strength of the resonator to the reservoir, such that Ylk ~^ 
(27r)^^ Jg°° dyucr (yu), we have to solve integrals appearing in Eq. ([3]), related to correlation 
functions of the form 



Jo 

g ^-^(,-.oKt-t') ^(^)Ar(^)A(/., t)A(/i, 0, (4) 

where the thermal average excitation of the reservoir N{fi) is defined by {b^fi)b{fi')) = 
N{fi)6{fi — fi'), while the system-reservoir coupling is modeled as 

with the parameter ^ accounting for the spectral sharpness around the TD frequency of the 
nonstationary mode. It is quite reasonable, for the case of weak system-reservoir coupling 
considered here, to assume a Lorentzian shape for the function A(/i, t), centered around the 
frequency uj{t). Moreover, as expected, an estimate of the time average of the operator A(t) 
reveals that the TD system- reservoir coupling falls with Xq/ \fi — ujq\, so that the larger the 
detuning, the smaller the coupling. Performing the variable transformations t = ( {t — t') 
and u = {ujq — fi) — sin(^t) in Eq. (jlj) and assuming, as usual, that a and are functions 
that vary slowly around the frequency ujq, we obtain 

Kt 

dt' {K\t)A{t')) = kk\N{ujq) / rfre-*"^(") 

Jo 

dv e^^^^ 

where, apart from the functions F{j) = cos {(t — t) + cos {(t) — r sin {(t), G{t) = sin {(t) — 
sin {(t — r), and a = ti;o/x~sin {(t), we have defined the dimensionless parameters x = Fo/C, 
K = ^/x, and s = xKi where Fq = cr(ct;o)AQ is the well-known damping rate of a stationary 
mode. Under the assumption that x/^o ^ the upper limit a can be extended to infinity 
and the corresponding integral can be evaluated analytically, leading to the correlation 
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function 



dt' {A^{t)A{t')) = 2N{ujQ)xn\ / dr 

Jo 



Ct ie[F(r) + irG(r) 



G'3(r)(l + e 



X <^ G{t) cos — + 2Ksin ' ^ 



= Nicuohit). (7) 

where = 2k/G{t) and 7(t) is related to an effective time- dependent damping rate. For 
the sake of completeness, before analyzing the influence of the parameters x, n, and e on 
the damping rate of a nonstationary mode, we compute its reduced density operator. To 
this end, assuming a reservoir at absolute zero, where N{ujq) = 0, we obtain from Eq. ([3]) 
the master equation 

= 2 Re [7(t)] ap{t)a'^ - -i* {t) p{t)a) a - -^{t)a^ap{t), (8) 

whose c-number version, for the normal ordered characteristic function xiViV*-:^) = 
Tr exp(?7a"'') exp(— ?7*a)] , is given by 

dt = ^'^"^^-^ dv* ■ 

Assuming a solution of the form xiVjV*^^) = x{'n{^)iV*{^))^ we obtain rfit) = ?7o e"'"'^*^/^, 
where rjo = rj{t = 0) and T(t) = J^'j{T)dT is the effective damping rate. Assuming, in 
addition, that xiVyV*^^) = xl'?; "'?*) ^ = 0)l,,^,,(t)) "we obtain from the Glauber-Sudarshan 
P-representation and the initial superposition state |\E'(t = 0)) = J2e^i\^oi), the reduced 
density operator of the nonstationary mode 

p(i)=Ar25^^^,C«,(t)|«.(t))(«Ht)l, (10) 
where a£{t) = ao£e~'^^^'' and 



Cu'it) = exp 



I {\ao,f + M') + a*,,aoi [l - e'^^^PWl] | 4q. (11) 

We note that, as expected, the decay rate turns out to be a real function even when r(t) is 
complex. For the particular case where the nonstationary mode is prepared in the superpo- 
sition state |\E'(0)) = Af {\ao) + |— ao)), the function multiplying the nondiagonal elements 
of the density matrix reads 

Ci2(t)=exp[-2|aor(l-e-2^'=([^Wl)]- (12) 
5 



We now analyze the influence of the parameters x, k, and e on the effective damping 
rate r{t) which, in its turn, determines the decoherence time of the superposition |\E'(0)), 
as given by Eqs. fill I) and (1121) . Starting with the parameter x = Fq/C, a measure of the 
rate of variation of the frequency compared to the damping constant Tq, it is evident 
from Eq. (j7]), as expected, that the damping function T(t) decreases in proportion to x. 
Otherwise, in the adiabatic regime where ( approaches Fq we also expect F(t) to be close to 
the damping constant Fq. Regarding parameter k = ^/x? which accounts for the range of 
oscillation of uj{t) compared to the Lorentzian sharpness ^, we expect the damping function 
to decrease as the range of oscillation x increases, as long as the variation rate ( is adjusted 
to be significantly higher than Fq. When both parameters k and x are adjusted accordingly, 
to be significantly smaller than unity, the system-reservoir coupling is weakened as well as 
the damping function F(t), consequently increasing the decoherence times of superposition 
states. Differently from x, our expectation concerning k, is blurred in Eq. ([7j): just as it is 
confirmed by the factor k,^, it is refuted by the decay function e~'^'^^ in the integral. Finally, 
the parameter e = x/C ^^^ly also be defined as e = x/k, as long as the damping rate Fq 
approximates the sharpness ^, weighting the contributions of parameters x and k. For the 
same reason as k, the role played by e in the behavior of F(t) is also blurred in Eq. ([7]). 

To clarify the role of the parameters x and k in the damping rate, in Fig. 1(a) we 
plot the function Ci2(t) against the scaled time Fgi, considering the initial superposition 
|\E'(0)) = Af {\ao) + \—ao)) with ao = 1. The thick sohd line corresponds to the case of a 
stationary mode where uj(t) = uq, prompting the expected result T(t) = Tot/2. Setting k = 
1/2, the solid and dashed lines correspond to x = 1/2 and 1/10, respectively. As expected, 
the damping function decreases as the rate of variation of the frequency increases. In fact, 
a higher rate of variation works to hinder the system-reservoir coupling, lengthening the 
response time of the system. With k = 1/10, the dashed-dotted and dotted lines correspond 
to X = 1/2 and x = 1/10, showing that the amplitude of oscillation x is more effective in 
diminishing the damping rate than the rate of variation (. This unexpected result reveals 
interesting aspects of the physics of nonstationary cavity modes: first of all, as demonstrated 
below, i) the time-dependence of uj(t) - the values of the frequencies x C ^ required 
to practically switch off the system-reservoir coupling can be engineered through atom-field 
interaction; furthermore, ii) in the adiabatic regime, where C/^o, x/^o 1; the atom-field 
interaction is still modelled by the Jaynes-Cummings interaction despite the nonstationary 
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mode [3]. Consequently, all the protocols developed for the implementation of processes 
in stationary modes — for example, quantum state or Hamiltonian engineering and logical 
devices — become directly applicable to the nonstationary mode considered here. 

Figs, l(b-f) display the damping process in the evolution of the amplitude of the coherent 
state a(t) = aoexp[—iQ(t) —T{t)] composing the superposition 1^1' (t)). All these figures 
refer to the same time interval as that used in Fig. 1(a), thus leading to the same number 
of cycles coming from the rotation in phase space, due to the factor e~*^(*). In all figures 
the ratio ooq/Tq = 10 is set to a fictitious scale to make clear the spiraling of a{t). In Fig. 
1(b), related to the thick solid line of Fig. 1(a), we observe the loss of excitation carrying the 
initial coherent state to the vacuum state. In this figure we also plot, in a thick dotted line, 
the evolution of the amplitude —a{t) of the other component of the superposition state. 
Figs, l(c-f) correspond respectively to the solid, dashed, dashed-dotted and dotted lines of 
Fig. 1(a), showing a gradual suppression of the loss of excitation which, differently from Fig. 
1(b), does not occur at a uniform rate, due to the oscillatory character - coming from Eq. 
([7]) - of their corresponding curves in Fig. 1(a). Fig. 1(e) clearly reveals this nonuniform 
character of the excitation loss through the distinct gaps between the cycles described by 
the amplitude a{t) on its (obstructed) way to the vacuum. 

Next, considering some sensitive features in the present scheme to control decoherence, 
we first address the time-dependent system-reservoir coupling Afc(t), which can be justified 
through the treatment of two coupled harmonic oscillators, one of them with time- dependent 
frequency. We start with the usual coupling CX1X2, where Xi(t) = Ci(t)(ai + a[) and 
X2 = C2(a2 + al). Within the interaction picture and the rotating- wave approximation, we 
end up with an time-dependent interaction of the form C(t) (a\a2 + aia\], similar to what 



had been considered in Refs. 



13 



19 



20I, 



2l|. Since a Lorentzian function applies whenever 



we have weak system-reservoir coupling, the time- dependent function assumed in Eq. ([5]) 
follows straightforwardly. 

The most sensitive point, however, is the engineering of the nonstationary mode whose 
state is to be protected. There is a great deal of literature exploring nonstationary modes. 



especially in respect of Casimir effect 



22] . We present below a scheme to engineer a nonsta- 



tionary mode uj{t) = cuq + x sin(C^) from the interaction of a driven two- level atom (frequency 
uja) with a stationary cavity mode (frequency oOc) given by 
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H = uj^a^a + y(T, + Fit) (a+e"^^^* + a_e^^^*) + G {aa+ + aV_) , (13) 

where ujl stands for the frequency of the classical driving field and G denotes the Rabi 
frequency. The atomic operators are given by o"^ = |e) (e| — \g) {g\, cr+ = |e) {g\, and 
^- = Ifl') (^l, e and g being the excited and the ground states. We assume the atomic ampli- 
fication modulated as F{t) = FQCOs{(t/2 + 0). In the interaction picture, the transformed 
Hamiltonian is given by 



Hi = G {aa+e'^'' + aV-c"'^^*) + F{t) {a+e'^'' + a.e-'^'') , (14) 

where 6i = Ua — oj^ and 82 = uja — ojl are the atom-field and the atom-laser detunings. 
Next, we define Hi = G {aa+e'^^^ + a'^a-e"'^'^) and H2 = F{t) {a+e'^^* + a^e~'^^^) , and 
assume the highly off-resonance laser amplification process, such that \62\ » FoX,G,\6i\ 
with G << \6i\. Under this assumption, the strongly oscillating terms of H2 lead, to a good 
approximation, to the effective Hamiltonian {2^, 

n,ff = Hi-tH2{t) j H2{T)dT 

= Uca^a + fi(t)o"2 + g (a(7+ + aV_) , (15) 
where Q{t) = uJa/2 + F'^{t)/Si. The diagonalization of Hamiltonian Heff is easily accom- 



phshed through the dressed atomic basis {\g,n),\e,n — 1)} [2J]. Under the usual assumption 
that G'^n -C 5^, we obtain the dispersive atom- field interaction: 

n = u^a^a + VL{t)a, + i:{t)a^aa,, (16) 

where the adjustment = 7r/4 makes T(t) = Ti + T2 sin((^t) with Ti = 
[1 - 3F2/25i52] GV^i and T2 = (CV^i) {F^/'25i52). Evidently, by turning off the laser we 
obtain the usual Stark shift Ta^aaz with T = G^/6i. In a frame rotating with the shifted 
atomic frequency fi(t), obtained through the unitary operator U = exp —iQ{t)az with 

n(t) = J n{t')dt', the state vector associated with the transformed Hamiltonian H = 
ujcO^a + T[t)a^aaz is given by 

|vi> {t)) = e^^W \g) 1$, [t]) + e-^W |e) |$e [t]) , (17) 
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where, in the Fock basis: (t)) = ^„ {C-,n\ \1/ {t)) |n), £ = g,e. Using the orthogonahty 
of the atomic states and Eqs. f[T6|) and f|T7|) . we obtain the uncoupled TD Schrodinger 
equations 



= co'^(t)a"'"a, 



(18) 
(19) 



where Ug = Uc — T(t) and Ue = uJc + 'Y{t)- Therefore, when preparing the atom in the 
fundamental state, we obtain the TD frequency u!{t) = cuq^X sin(Ci) where, from interaction 
(fT6|) . ujq = ujc + Ti and x = ^2- Note that the atom crosses the cavity remaining in its 
ground state (due to its off-resonance interactions with both the cavity mode and the classical 
field), so that there is no injection of noise coming from the atomic decay to the engineered 
nonstationary cavity mode. Assuming typical values for the parameters involved in cavity 



QED experiments 



25 



26]: G ~ 3 X 10V\ \5i\ ~ 10V\ l^sl ~ 10V\ and Tq ~ lOVS 
it follows, with the intensity Fq ~ 10 x G, that x ~ 4 x lO^s^^ with ( < lO^s"^. Since it is 
reasonable to assume ^ ~ To, the value 1/10 for the parameters k and x employed to obtain 
the dotted line of Fig. 1(a), is easily accomplished. 

Evidently, to circumvent the difficulties introduced by the small time interval of atom- 
field interaction, it would be interesting to engineer the nonstationary mode through a 
sequential interaction of atoms, one by one, with the cavity mode. The trapping of an atom 



inside the cavity, along the lines suggested in Ref. [27|, is also a possibility to be analyzed. 
Otherwise, nonstationary modes can also be achieved by other schemes as the mechanical 
motion of the cavity walls {28 1, suitable for our purpose since the frequency attainable is in 
the gigahertz range, or even more sophisticated schemes where the effective motion of the 
walls is generated by the excitation of a plasma in a semiconductor 29l|." 

We have thus presented in this paper a scheme which practically switches off the reservoir 
of a cavity field by engineering a suitable nonstationary mode uj{t). Besides analyzing 
the physical parameters required to accomplish this process, we also demonstrated how to 
engineer such a nonstationary mode through its dispersive interaction with a driven atomic 
system. Evidently, the scheme presented here for a time- dependent cavity mode applies 
directly to any oscillatory system such as trapped ions, nanomechanical oscillators, and 
superconducting transmission lines; it can also be extended to any nonstationary quantum 
system. We believe that both techniques presented here, to protect quantum states through 



nonstationary quantum systems and to engineer such systems, can play an essential role in 
quantum information theory. 

Beyond the information theory, we believe that the present work can directly contribute 
to the field of cavity quantum electrodynamics, specifically in the less-explored topic of the 
interaction of a two-level atom with a nonstationary mode in the adiabatic regime. In fact, 
as the engineering of a nonstationary mode is relatively easy to be accomplished in the 



adiabatic regime — through the mechanical motion of the cavity walls [28| or atom-field 
interaction, as demonstrated in this work — typical quantum optical phenomena may be 
investigated in this particular context. 
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Fig. 1. In (a) we plot the function Ci2{t) against the scaled time Tot. Considering a plot 
of lm{a{t)) against Re(Q;(t)), in (b-f) we observe the damping process in the evolution of 
the amplitude of the coherent state a{t). 
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